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Abstract 

The spin-one-half FaHcov-Kimball model with spin-dependent on-site in- 
teraction between localized (/) and itinerant (d) electrons is studied by small- 
cluster exact-diagonalization calculations and a well-controlled approximative 
method in two dimensions. The results obtained are used to categorize the 
ground-state configurations according to common features (charge and spin 
ordering) for all / and d electron concentrations (nj and n^) on finite square 
lattices. It is shown that only a few configuration types form the basic struc- 
ture of the charge phase diagram in the nj — rid plane. In particular, the 
largest regions of stability correspond to the phase segregated configurations, 
the axial striped configurations and configurations that can be considered as 
mixtures of chessboard configurations and the full (empty) lattice. Since the 
magnetic phase diagram is much richer than the charge phase diagram, the 
magnetic superstructures are examined only at selected values of / and d elec- 
tron concentrations. 

PACS numbers.:75.10.Lp, 71.27.+a, 71.28.+d 
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1 Introduction 



The interplay between charge and spin degrees of freedom in strongly correlated sys- 
tems has triggered enormous interest in recent years due to the rich variety of charge 
and spin orderings found in some rare-earth and transition-metal compounds. Charge 
and spin superstructures have been observed, for example, in doped niclate [1], 
cuprate [2] and cobaltate [3] materials, some of which constitute materials that ex- 
hibit high-temperature superconductivity. One of the simplest models suitable to 
describe charge-ordered phases in interacting electron systems is the Falicov-Kimball 
model (FKM) Indeed, it was shown that the simplest version of this model (the 
spinless FKM) already exhibits an extremely rich spectrum of charge ordered solu- 
tions, including various types of periodic, phase-separated and striped phases [5], [6]. 
However, the spinless version of the FKM, although nontrivial, is not able to ac- 
count for all aspects of real experiments. For example, many experiments show that 
a charge superstructure is accompanied by a magnetic superstructure [U El E] • In 
order to describe both types of ordering in the unified picture, a simple model based 
on a generalization of the spin-one-half FKM with an anisotropic, spin-dependent in- 
teraction that couples the localized and itinerant subsystems was proposed ^ . Thus 
the model Hamiltonian can be written as 

H =Y1 ^iAdja + f^fiadf^,di^> + J "^{f^^fi-a " f ^ fia) df^di^ , (1) 

ijcr iaa' icy 

where /J^, are the creation and annihilation operators for an electron of spin 
a =tj i in the localized state at lattice site i and c?^, di^j are the creation and anni- 
hilation operators of the itinerant electrons in the d-band Wannier state at site i. 

The first term of (1) is the kinetic energy corresponding to quantum- mechanical 
hopping of the itinerant d electrons between sites % and j. These intersite hopping 
transitions are described by the matrix elements tjj, which are —t if i and j are the 
nearest neighbours and zero otherwise (in the following all parameters are measured 
in units of t\ The second term represents the on-site Coulomb interaction between 
the (i-band electrons with density rid, = ^dl^ = \ Yl,ia dtadia and the localized / 
electrons with density n/ = Nj /L = J2ia fiafi^i where L is the number of lattice 
sites. The third term is the above mentioned anisotropic, spin-dependent local in- 
teraction of the Ising type between the localized and itinerant electrons that reflects 
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the Hund's rule force. Moreover, it is assumed that the on-site Coulomb interaction 
between / electrons is infinite and so the double occupancy of / orbitals is forbidden. 

Since the /-electron occupation number f^fia of each site i still commutes with 
the Hamiltonian (1), the /-electron occupation number is a good quantum number, 
taking only two values: Wi^j = 1 or 0, according to whether or not the site i is 
occupied by the localized / electron. Now the Hamiltonian (1) can be written as 

h=y: K.did,^, (2) 

ijcr 

where hij = Uj + {Uwi + Jwi-a — Jwia)5ij. Thus for a given /-electron configuration 
w = {wi,W2, . . . ,wl}, the Hamiltonian (2) is the second-quantized version of the 
single-particle Hamiltonian h{w), so the investigation of the model (2) is reduced to 
the investigation of the spectrum of h for different configurations of / electrons. This 
can be performed exactly, over the full set of /-electron distributions (including their 
spins), or approximatively, over the reduced set of /-electron configurations. Here we 
use a combination of both numerical methods. The same procedure has been used 
already in our previous paper [9j to study the ground-states of the extended FKM 
in one dimension. In the mentioned paper we have also presented some preliminary 
results concerning the ground states of the model in two dimensions. These results 
revealed that the extended FKM with spin-dependent interaction between / and d 
electrons can describe various types of charge and magnetic superstructures. In the 
present paper we try to construct the comprehensive phase diagram (in the Nf — 
plane) of the extended FKM in two dimensions. We supply our studies by a detailed 
finite-size scaling analyses in order to minimize the infiuence of finite-size effects on 
the ground-state properties of the model. Here we consider only the case of strong 
Coulomb interactions {U = 4), since the one-dimensional studies showed that the 
influence of the Hubbard interaction (between d electrons) on ground-states of the 
model Hamiltonian (2) can be neglected in this limit. At the end of this paper we 
specify precisely conditions under which this term can be neglected in two dimensions. 

2 Results and discussion 

To construct the comprehensive picture of charge and magnetic ordering in the ex- 
tended FKM in two dimensions, the complete phase diagram of the model in the 
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Nf — Nd plane has been calculated point by point for all even number of Nf and 
Nd. Of course, such a procedure demands a considerable amount of CPU time, that 
imposes severe restrictions on the size of clusters that can be studied numerically 
(L = 8 X 8). First we have concerned our attention on the problem of charge order- 
ing. In Fig. 1 we present results of our numerical calculations obtained for U = 4 
and J = 0.5 in the form of the skeleton phase diagram. One of the most interesting 
observations is that the phase diagram consists of only a few configuration types, 
although the total number of possible configurations increases very rapidly with the 
cluster size L as 3^. In particular, we have detected 5 different charge configuration 
types, and namely: (a) the segregated configurations, where / electrons clump to- 
gether, (b) the n-molecular phases, which have been observed only for small N^, and 
small Nf, (c) the axial stripes, (d) the regular phases (stable only in isolated points) 
and mixtures of regular (usually chessboard) phases and empty/full lattice accom- 
panied by (e) the miscellaneous configurations. The typical examples corresponding 
to these configuration types arc depicted in the lower part of Fig. 1. As one can see 
the largest stability regions correspond only three configuration types, and namely, 
the segregated configurations and the axial stripes, which fill the left and right side 
of the phase diagram and the regular phases/mixtures of regular phases and empty 
(full) lattice locahzed at the central part of the skeleton phase diagram. Moreover, it 
was found that the homogeneous /-electron distributions (/ electrons are distributed 
as far away from each other as possible) are the ground states only in the region d 
and only in isolated points, while in the rest part of the phase diagram the inho- 
mogeneous distributions (the segregated phases, the axial stripes, the n-molecular 
phases) are preferred. This result is very valuable since in the past years a strong 
interest in the experimental and theoretical studies of strongly correlated systems 
was focused on the physics that leads to an inhomogeneous ordering, especially to an 
inhomogeneous charge stripe order due to the observation of such ordering in doped 
niclate and cuprate materials some of which exhibit high-temperature superconduc- 
tivity. From the theoretical point of view, the presented skeleton phase diagram 
clearly demonstrates that this relatively simple model (generalized two-dimensional 
FKM) can describe such inhomogeneous stripe ordering. Moreover, it was shown 
that the stability area of stripes is relatively large and, in addition, includes both 
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physically interesting conditions Nf + — L and Nf + — 2L. 

Let us now turn our attention on the second problem, and namely, the prob- 
lem of spin ordering. We have found, that although the charge phase diagram of 
the generalized two-dimensional FKM is rather simple, the spectrum of magnetic 
solutions is very rich. As demonstrated in Fig. 2, the one charge distribution could 
have many different spin configurations and therefore the exact classification is very 
difficult. Moreover, in the two-dimensional case the finite size effects on spin order- 
ings are still large for clusters with L < 64, and thus we were not able to construct 
definitive picture of magnetic phase diagram for all Nf and N^. For this reason we 
have focused our attention on the physically the most interesting cases mentioned 
above, and namely, Nf + Nd — L and Nf + Nd — 2L, where exhaustive studies have 
been performed. To minimize the finite size effects we have studied the set of finite 
clusters of L = 6 x 6, 8 x 8, 10 x 10, 12 x 12, 16 x 16, 18 x 18 and 20 x 20 sites. 

We have started our study with the case Nf + Nd = 2L, which is slightly simple 
for a description. As shown in Fig. 3 the ground states for Nf + A^^^ = 2L are 
antiferromagnetic (AF) with alternating pattern, where the electrons (for n/ < 1/2) 
or holes (for Uf > 1/2) form the axial distributions. Our calculations showed that 
these inhomogeneous stripe distributions are stable for large Coulomb interactions, 
while decreasing U leads to their destruction and prefers the homogeneous electron 
arrangement (see Fig. 4). 

From the skeleton phase diagram we know that the condition line Nf + Nd — L 
lies in the axial stripe area, but in comparison with the previous case the situation 
is fully different. The first fundamental difference is that for sufficiently small /- 
electron concentrations Uf the ground state could be ferromagnetically (F) ordered 
(sec Fig. 5). The second one is, that although with increasing /-electron concentra- 
tion the ground states are the AF, these AF arrangements are formed by F ordered 
clusters (domains). In addition, for Uf — 1/4 and Uf — 1/3 (L > 144) a new 
type of stripes (known as the ladders) is occurred. And finally, a detailed analyses 
showed that there exists a critical /-electron concentration ~ 1/4 bellow which 
the ground-states are phase separated. 

To verify the ability of our method to describe ground-state properties of macro- 
scopic systems we have performed similar calculations for U — 8 and several values 
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of Uf and {nj = 3/4, Ud = 1/2; n/ = 2/3, = 2/3; Uf = 1/2, = 1; n/ = 1/3, 
nrf = 4/3; = 1/4, rid = 3/2) for which there exist the numerical results obtained in 
the thermodynamic limit {L oo) using the method of restricted phase diagrams [8]. 
The numerical results for ground-state configurations obtained by these two differ- 
ent approaches are compared in Fig. 6. It is seen that both methods yield the same 
results for the charge distributions of /-electrons, but differ in a prediction of spin 
distributions. To check the stability of our solutions we have performed numerical 
calculations also for L = 12 x 12, 16 x 16 and 18 x 18 clusters. The results obtained are 
shown in Fig. 7 and they clearly demonstrate that ground states found for L = 8 x 8 
cluster hold also on much higher clusters. This fact allows us to extrapolate our 
numerical results to the thermodynamic limit, where they can be directly compared 
with the Lemanski's results. This comparison shows (see Fig. 8) that our method 
yields in all cases the lower energy in the thermodynamic limit {L oo) than the 
Lemanski's approach. Moreover, a comparison of our and Lemanski's ground states 
(Fig. 6) provides a direct explanation of this discrepancy, and namely, that the unit 
cells used in Ref. jS] are too small to describe correctly the spin distributions. 

Although we have presented here only the basic types of charge and magnetic 
superstructures they clearly demonstrate an ability of the model to describe dif- 
ferent types of charge and magnetic ordering. This opens an alternative route for 
understanding of formation an inhomogeneous charge/ magnetic order in strongly 
correlated electron systems. In comparison to previous studies of this phenomenon 
based on the Hubbard [lO] and t — J model [11], the study within the generalized 
spin-one-half FKM has one essential advantage and namely that it can be performed 
in a controllable way (due to the condition [f^fia-,H] = 0), and in addition it al- 
lows easily to incorporate and examine effects of various factors (e.g., an external 
magnetic field, nonlocal interactions, etc.) on formation of charge and magnetic 
superstructures. 

Of course, one can ask if these results persist also in the more realistic situations 
when additional interaction terms are included into the Hamiltonian (2). From 
the major interaction terms that come into account for the interacting d and / 
electron subsystems only the Hubbard type interaction Udd J2i df^di^df^di^ between 
the spin-up and spin-down d electrons has been omitted in the Hamiltonian (2). 
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In work [8] Lemanski presents a simple justification for tlie omission of this term, 
based on an intuitive argument: the longer time electrons occupy the same site, 
the more important becomes interaction between them. According to this rule the 
interaction between the itinerant d electrons {Udd) is smaller than the interaction 
between the localized / electrons (f///) as well as smaller than the spin-independent 
interaction between the localized and itinerant electrons. Here we specify more 
precisely conditions when this term can be neglected. To determine the effects of 
Udd interaction on the ground-states of the spin-one-half FKM with spin-dependent 
interaction (J = 0.5) in two-dimensions the exhaustive studies of the ground-state 
phase diagram of the model (in the nj — Udd plane) have been performed. Of course, 
an inclusion of the Udd term makes the Hamiltonian (2) intractable by methods used 
for the conventional spin-one-half /spinless FKM and thus it was necessary to use 
other numerical methods. Here we used the Lanczos method [12] to study exactly 
the ground states of the spin-one-half FKM generalized with Udd interaction between 
the spin-up and spin-down d electrons. Such a procedure demands in practice a 
considerable amount of CPU time, which imposes severe restrictions on the size of 
clusters that can be studied within the exact-diagonalization method. For this reason 
we were able to investigate exactly only the clusters up to L = 10. The results of 
numerical calculations obtained for f/ = 4 are summarized in Fig. 9 in the form of 
Uf — Udd phase diagram (the half- filled band case nj + nd = 1 is considered). One 
can see that the ground-state configuration w^{Nf) found for Udd = persists as a 
ground state up to relatively large values of Udd iU^j. ~ 3.5), revealing small effects 
of the Udd term on the ground states of the model in the strong U interaction limit. 
For this reason our numerical calculations have been done exclusively for large U . 

In conclusion, we have used a combination of small-cluster exact diagonaliza- 
tion calculations and a well-controlled numerical method to study the ground-state 
properties of the spin-one-half FKM extended by spin-dependent on-site interaction 
between localized and itinerant electrons in two dimensions. The results obtained 
have been used to construct the comprehensive picture of charge and spin ordering 
in this model. It was shown that only a few configuration types form the basic 
structure of the charge phase diagram in the Nf — Nd plane. In particular, the 
largest regions of stability correspond to the phase segregated configurations, the 
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axial striped configurations and configurations that can be considered as mixtures of 
chessboard configurations and the full (empty) lattice. Moreover, it was found that 
the model exhibits a rich spectrum of magnetic solutions including various types of 
ferro- and antiferromagnetically ordered phases. 

This work was supported by Slovak Grant Agency VEGA under Grant No.2/7057/27 
and Slovak Research and Development Agency (APVV) under Grant LPP-0047-06. 
H.C. acknowledges support of Stefan Schwartz Foundation. 
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Figure 1: The skeleton phase diagram of the two-dimensional spin-one-half FKM 
extended by spin-dependent interaction calculated for L = 8 x 8, C/ — A and J = 0.5. 

(a) The segregated configurations, (b) the n-molecular phases, (c) the axial stripes, 
(d) the regular phases, the mixtures of regular phases and full/empty lattice and (e) 
the miscellaneous phases. Lower part: the typical examples of ground states of the 
model. Large dots: sites occupied by / electrons, small dots: vacant sites. 



10 



N=32,N =16 
f d 



N=32,N^=16 



N=32,N =110 
f d 



N=32,N =112 
f d 



AAAAAAAA AA^^AAy^ ▼▼▼AyAAA 

AAAAAAAA y^AAYYAA ▼AyA^AvA 

TTTTTTTT ^^TT^^TT TAyAvA^A 

TTTTTTVV TY^^YY^^ TTT^T^^^ 

N=32,N=114 N=32,N=116 N=32,N=120 

f ' d f d f d 



N=36,N =62 N=36,N =62 N=36,N =64 

f d f d f d 

•••••••• -T-T-T-T .A-v-A-A 

• ▼•▼•A-A. ▼•▼•A.^. 
•••••••• •▼▼A-T-A .AAA.A.^ 

• ••••••• yr ' ^ Y Y ' T ' A-yAA-y. 

• ••••••• ' A ' A A Y ' ^ •Y-yAy-A 

•••••••• y>A>yAA> Y>y<AAA- 

• ••••••• •Y>Y>A>A •A>A<Y<A 

• ••••••• A>A>A>A> Y<Y<Y<Y> 



▼AvAvAyA yyyyyyyy AA^vAAfv 

▼AyAyAyA AAAAAAAA AAyT^AyT 

YAyAyAyA yttttttt ttaaytaa 

YAyAyAyA AAAAAAAA ytAAytAA 

N=32,N =40 N=32,N =40 N=32,N =88 N=32,N =96 

f d f d f d f d 



YYTTAAAA AAyAytAy TAyAyAyA 

AAAAyTTT TTAyAAyA AyAyAyAy 

YYTTAAAA AAyAytAy tAyAyAyA 

AAAAyTTT TTAyAAyA AyAyAyAy 



Figure 2: Several spin distributions with common charge pattern depicted on L = 
8x8 cluster. To visualize spin distributions we use A for the up spin electrons and 
V for the down spin electrons. 
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Figure 3: Typical ground-state configurations of the two-dimensional generalized 
FKM obtained for selected values of n/ on finite clusters of L = 8 x 8, L = 10 x 10 
and L = 12 X 12 sites at C/ = 4, J = 0.5 and n/ + = 2. 
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Figure 4: Ground-state configurations of tlie two-dimensional generalized FKM ob- 
tained for different U {U = A and U = 2) and Uf {uf = 1/4, Uf = 1/8 at Uf + na = 2). 
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Figure 5: Typical ground-state configurations of the two-dimensional generalized 
FKM obtained for selected values of on finite clusters of L = 8 x 8, L = 10 x 10 
and L = 12 X 12 sites at C/ = 4, J = 0.5 and Uf + nd^ 1. 
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Figure 6: The ground-state configurations of tlie two-dimensional spin-1/2 FKM for 
[7 = 8 and four different pairs oi nj and obtained by two different approaclies: 
tlie metliod of restricted pliase diagrams ^ (-D1, D2, D4 and Db) and our numerical 
metliod {D1',D2',DA' and D5'). Tfie sliaded region in tlie lower left corner shows 
the unit cell, and line segments show the translation vectors that are used to tile the 
two-dimensional plane. 
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Figure 7: The ground-state configurations of the two-dimensional spin-1/2 FKM 
calculated for the same U^Uf and rid values as in Fig. 6 on different clusters of 
L = 8 X 8, 12 X 12 and 16 x 16 sites. 
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Figure 8: The difference between our and the Lemanski's ground-state energies as 
a function of 1/L calculated for the extrapolated configurations corresponding to 
Dl', D2', DA', D5' and Dl, D2, DA, D5 phases. The fines are only guides to the eye. 
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Figure 9: The ground-state phase diagram of the spin-one-half FKM extended by 
the Hubbard interaction between the itinerant electrons calculated for [/ = 4 and 
J = 0.5 on small finite cluster of L = 10 sites. Below U^^ the ground states are the 
ground-state configurations of the conventional spin-one-half FKM {Um = 0). Above 
these ground states become unstable. The two-dimensional exact diagonalization 
results. 
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